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' The parametric correlations of the transmission eigenvalues Tj of a A^-channel 

quantum scatterer are calculated assuming two different Brownian motion 
^ . ensembles. The first one is the original ensemble introduced by Dyson and 

assumes an isotropic diffusion for the S'-matrix. We derive the correspond- 
I ing Fokker-Planck equation for the transmission eigenvalues, which can be 

' mapped for the unitary case onto an exactly solvable problem of N non- 

interacting fermions in one dimension with imaginary time. We recover for 



o 



0^ . the Tj the same universal parametric correlation than the ones recently ob- 

^ , tained for the energy levels, within certain limits. As an application, we con- 

g I sider transmission through two chaotic cavities weakly coupled by a n-channel 

I ' point contact when a magnetic field is applied. The 5-matrix of each chaotic 

■ cavity is assumed to belong to the Dyson circular unitary ensemble (CUE) 
O ■ and one has a 2x CUE — > one CUE crossover when n increases. We calculate 
^ ' all types of correlation functions for the transmission eigenvalues Tj and we 

get exact finite results for the averaged conductance {g) and its variance 
^ ' {6g'^), as a function of the parameter n. The second Brownian motion en- 

■ semble assumes for the transfer matrix M an isotropic diffusion yielded by a 
multiplicative combination law. This model is known to describe a disordered 
wire of length L and gives another Fokker-Planck equation which describes 
the L-dependence of the Tj. An exact solution of this equation in the unitary 
case has recently been obtained by Beenakker and Rejaei, which gives their 
L-dependent joint probability distribution. Using this result, we show how to 
calculate all types of correlation functions, for arbitrary L and N. This allows 
us to get an integral expression for the average conductance which coincides 
in the limit N ^ oo with the microscopic non linear c-model results obtained 
by Zirnbauer et al, establishing the equivalence of the two approaches. We 
review the qualitative differences between transmission through two weakly 
coupled quantum dots and through a disordered line and we discuss the math- 
ematical analogies between the Fokker-Planck equations of the two Brownian 
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motion models. 
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I. INTRODUCTION 

It is now rather well established that the quantum energy levels of classically chaotic 
billiards have Wigner-Dyson correlations ||l],3. These universal correlations are closely as- 
sociated with the concept of "quantum chaos" and have first been obtained in random matrix 
theory (RMT) in a different context (statistical description of the spectral properties of 
complex nuclei and small metallic particles). Futhermore, many works p-pD| have consid- 
ered the case where the system Hamiltonian H depends on an external parameter t. Instead 
of a single matrix ensemble defined by a certain probability density, a continuous family of 
ensembles, defined by a t-dependent probability density, is considered. The universality of 
the level correlations can then be extended to a broader domain, which involves not only 
(small) level separations at a given t, but also (small) external parameter separation 5t. This 
broader universality has been numerically checked in very different physical systems: e. g. 
disordered conductors 0, quantum billiards |Tl|], the hydrogen atom in a magnetic field |1^ 
and correlated electron systems |ll3| , p!^ . Analytical derivations use perturbation theory 
or the non linear a-model [|^. Another fruitful way to calculate the universal expressions of 
those parametric level correlations is based on Brownian motion ensembles of random ma- 
trices, first introduced by Dyson. The idea is to assume that the N x N Hamiltonian matrix 
diffuses when one varies the parameter t in the available manifold given its symmetries. The 
t-dependence of the level distribution is then determined by a Fokker-Planck equation , 
which can be solved in the large A^-limit for arbitrary symmetries and for finite N in the 
unitary case (Sutherland method). 
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For quantum transport in the mesoscopic regime, the system is not closed, but open 
to electron reservoirs, and becomes a A^-channel scatterer with transmission eigenvalues Tj. 
Many transport properties |T5[ (conductance, quantum shot noise...) are linear statistics 
of the N eigenvalues Tj, which are the appropriate levels for a scattering problem. The 
scattering matrix 5* or transfer matrix M, suitably parametrized by the transmission 
eigenvalues Tj and certain auxiliary unitary matrices, are the relevant matrices for which 
RMT approaches have been formulated. For instance, if the unitary S'-matrix is distributed 



according to one of Dyson's classical circular ensembles ||22|, the corresponding transmission 
eigenvalue distributions have been recently obtained [|T3,rB|, exhibiting the same logaritmic 
pairwise repulsion than for the energy levels, and hence giving essentially the same Wigner- 
Dyson correlations for the transmission eigenvalues. This RMT description is mainly relevant 
for quantum transport through ballistic chaotic cavities, which is a recent field of significant 
experimental investigations. Such cavities can be made with semiconductor nanostructures 
known as quantum dots [T^ with a few channel contacts to two electron reservoirs. The 



validity of this RMT description is confirmed by microscopic semiclassical approaches |18 



and by numerical quantum calculations |T6|,[T9|. 

The universality of the transmission eigenvalue correlations are established for a given t, 
and therefore we want to know if it can be extended to their parametric dependence, too, 
as for the energy levels. Futhermore, we want to see if these parametric correlations obey a 
similar behavior than that of the energy levels. This is one of the issues which we address 
in this work, using two Brownian motion ensembles for the transmission eigenvalues, which 
result from two different assumptions concerning the t-dependence of S (S'-Brownian motion 
ensemble) or of M (M-Brownian motion ensemble). This issue has been recently considered 
also in Refs. 

The S'-Brownian motion ensemble is simply the original ensemble introduced by Dyson 
for the scattering matrix S 



s= ':,] ■ (1-1) 




The idea is to assume that 5* diffuses in the S'-matrix space with respect to a fictitious time 
t of the Brownian motion, and converges in the limit t — >■ cx3 to a stationary probability 
distribution given by one of the well known circular ensembles: circular unitary ensemble 
(CUE, (3 = 2) when there is no time reversal symmetry (applied magnetic field); circular or- 
thogonal ensemble (COE, (3 = 1) when there are time reversal and spin rotation symmetries 
( no spin-orbit coupling or applied magnetic field) and circular symplectic ensemble (CSE, 
/3 = 4) otherwise. Choosing an initial probability distribution at t = of higher symmetry 
than the stationary limit, e.g. two independent circular ensembles diffusing towards a single 
one, one can perform with the parameter t the crossover from the initial ensemble to the 
stationary ensemble. Actually, the Brownian motion ensembles COE CUE, CSE CUE, 
and 2xCUE — > CUE have been investigated in Ref. P^ , p5| in terms of the scattering phase 



shifts 9i. Semiclassical and numerical justifications of these models are given in Ref. p5[ 
The decisive difference between our approach and Refs. [0,^ is that we do not consider the 
usual eigenvector-eigenvalue parametrization of S, but the parametrization which explicitly 
uses the transmission eigenvalues Tj of tt'^ as introduced in Ref. [^. Since there exists 
neither a simple mathematical method nor an intuitive way to relate the 2N eigenvalues 
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e*^' of S with the transmission eigenvalues T,, it is justified to reconsider the Brownian 
motion in terms of the latter. This parametrization is suitable to determine the "time" 
dependence of the Tj, by a Fokker-Planck equation which we map, using a transformation 
introduced by Sutherland, onto a Schrodinger equation with imaginary time. In the case 
of the stationary CUE limit, one has a problem of non interacting fermions which we 
solve exactly for arbitrary A^. One obtains for the corresponding parametric correlations 
of the transmission eigenvalues the same universal form than for the energy levels, after an 
appropriate rescaling and within certain limits. 

A physical application of such a parametric ensemble can for example be realized by two 
weakly coupled ballistic (irregular) cavities, each of them connected to one electron reservoir 
(Fig. 1). The crossover between two CUE towards a single CUE is performed by changing 
the strength of the coupling from nearly uncoupled cavities (two CUE) to idealy coupled 
cavities (one CUE), when a small constant magnetic field is applied. The transmission 
eigenvalue parametrization is very well adapted to this case since the initial condition is 
easily realized by the requirement Tj = 0. Another application of the Dyson Brownian 
motion ensemble, concerning the crossover COE —>■ CUE in terms of the Tj is presented 



elsewhere 26 



Outside the problem of two weakly coupled quantum dots, we are also motivated by 



numerical calculations showing that the scattering matrix of a quasi one-dimensional 
disordered wire exhibits also some properties of the circular ensembles, in the sense that the 
correlation of the scattering phase shifts 6i were found to be approximately described by the 
universal correlations functions of these ensembles. For the disordered wire, a crossover from 
one circular ensemble towards two independent circular ensembles was observed as the length 



of the wire exceeds its localization length. The appropriate statistical description [p7H3T 



of the transmission eigenvalues of disordered wires is given by a different Brownian motion 
ensemble than the original Dyson model, and results from the multiplicative combination 
law of the transfer matrix M. This is what we call the M-Brownian motion ensemble 
which yields another Fokker-Planck equation taking into account the quasi one-dimensional 
structure. This Fokker-Planck equation has been solved by Beenakker and Rejaei in 
the unitary case, using the Sutherland transformation. A remarkable result which they 
found is that the transmission eigenvalue pairwise interaction is universal, in the sense that 
it does not contain any adjustable parameter. It coincides with the standard logaritmic 
RMT repulsion for small eigenvalue separations, and it is halved for larger separations. 
The physical interpretation of this halving is unclear, but must be somewhat related to the 
statistical decoupling of the reflection properties of the two opposite system edges. If it 
is right, the Dyson Brownian motion ensemble must also exhibit in the crossover regime 
2xCUE —>■ IxCUE a somewhat similar pairwise interaction for weak transmission. It is 
therefore an interesting idea to compare those two different Brownian motion ensembles for 
the transmission eigenvalues. Let us mention the obvious difference between them: their 
initial and stationary limits are inverse of each other. For the ^-matrix model, we start 
from no transmission towards a good CUE transmission. For the M-matrix model on the 
contrary, the initial condition corresponds to perfect transmission, and the system evolves 
with the time t (here to be identifled with conductor length L) to the zero transmission 
localized limit. 

This paper is composed of four main sections. First, we introduce the transmission 
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eigenvalue parametrization of 5* and define the S'-Brownian motion ensemble (section P A| ). 
The corresponding Fokker-Planck equation is derived in Appendix A. In section 



[TBi the 

Fokker-Planck equation is solved for an arbitrary initial condition in terms of fermionic one 
particle Green's functions, using the method introduced by Beenakker and Rejaei for the 

In section 111 CI, 



M-Brownian motion ensemble. 



we show how to recover for the transmis- 



sion eigenvalues the same universal parametric correlations than for the energy levels, after 
rescaling and within certain limits. 

In section |ITTt we treat completely the system of two coupled ballistic cavities by the 
S'-Brownian motion ensemble for the unitary case. The time of the Brownian motion is 
given in terms of the numbers of channels in the contacts to the electron reservoirs (A^) 
and in the contact between the two cavities (n). In Section [III A| , the crossover 2xCUE — >■ 
CUE is considered. The method of orthogonal polynomials, well known in the framework 
of random matrix theory 0, is extended to a more general crossover situation and yields 
the correlation functions as determinants (Eq. ( ^.18| )) containing a function Ki^{x,y]t) 
given as a sum of Legendre polynomials (Eq. (p.21|) ). This Section contains with Eqs. 
( p.27| - |3.30| ) also exact results for the average and the fluctuations of the conductance. In 
Section |III B| , we show that the transmission eigenvalue interaction deviates from the pure 
RMT logarithmic interaction, in the limit of weakly coupled cavities. One can note some 
interesting similarities with the pairwise interaction of the M-Brownian motion model given 
in Ref. 
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The third part of this paper (Section reconsiders the M-Brownian motion ensemble 
of Ref. ||2^ for disordered wires in the unitary case. We use the exact expression of the 
joint probability distribution found in Ref. and calculate the corresponding correlation 
functions, using a method very similar to the one used in Section III A . This is possible 
because, from a mathematical point of view, the Brownian motion models for S and M 
are very similar. We get an expression of the m-point correlations as a determinant (Eq. 
( [4.2(j| )) with a function K]si{X, A; t) given as a sum and an integral with Legendre polynomials 
and Legendre functions (Eq. ( [4.21| )). This allows us to prove the equivalence between the 
M-Brownian motion ensemble and a more microscopic approach, based on a non linear 
(T-model formulation and supersymmetry, as far as the behavior of the first two moments of 
the conductance (g) and {g'^) is concerned (section [IV B|) . We emphasize that the solution 
of the Brownian motion ensemble is now complete, contrary to the sigma model approach 
where one only knows (g) and {g"^) in the large A^-limit. In Section [IV Q , some further 
simplifications valid in the localized limit lead to asymptotic expressions for the density and 
the two point function for the logarithm of the transmission eigenvalues. 

In Section |V], we review some essential differences between the transmission properties 
characterizing two weakly coupled ballistic chaotic cavities (S-Brownian motion ensemble) 
on one side and an homogeneous disordered wire (M-Brownian motion ensemble) on the 
other side. We give some concluding remarks in section |V^, notably on the mathematical 
similarities between the S and M Brownian motion ensembles. 



II. 5 BROWNIAN MOTION ENSEMBLE FOR TRANSMISSION EIGENVALUES 



6 



A. Parametrization, definition and Fokker Planck equation 



The scattering matrix S* of a system connected to two electron reservoirs by two A^- 
channel contacts can be described in the transmission eigenvalue parametrization by 

/^i oWv/T^ iVT \(v2 o\ 

where vi^2i ^1,2 are N x N unitary matrices and T is a diagonal x A^-matrix with real 
entries < Tj < 1. The transmission matrix is t, but we mean by "transmission eigenvalues" 
Tj the eigenvalues of tH. This parametrization has recently ||l5| , p!6| been used to calculate the 
conductance properties of a ballistic cavity, which is known to be suitably described by one of 
Dyson's circular ensembles In the unitary case {f3 = 2) all unitary matrices vi^2, ^1,2 are 



independent whereas for the orthogonal {[3 = 1) and symplectic (/? = 4) symmetry classes, 
one has V2 = and U2 = (with = M'^ in the orthogonal case and = JM'^J'^, 
J = —i(Jy, in the symplectic case). In the following, the transmission eigenvalues Tj are 
described by angles ipj G [0, 7r/2] via Tj = sin^ ipj and therefore 





-P T (2.2) 



/ COS (p tsmip \ 
[isirnp cos (p j 

where is a diagonal matrix with entries (pj. The invariant measure for S in the parametriza- 
tion ( p.l|) and the change of variables (|2.2|) is given by ||T5| , [T6[| 



fi{dS) = F/^{{p)Y\_dipj fi{dui) fi{du2) fi{dvi) fj,{dv2) (2.3) 
j 

with 

Ff^iif) = n jsin^ "Pj - sin^ V^^f H (sm^"^ cos cpj) , /3 = 1, 2, 4 . (2.4) 

j<k j 

In the cases /? = 1,4 only two products ^{dui) fi{dvi) appear. 

The Brownian motion ensemble which we consider (cp. Dyson |^) describes a unitary 
random matrix which depends on a fictitious time. A small change on the time variable 
t ^ t + 6t yields (in the unitary case) for 5* the change: 

S{t + 6t) = S{t)e'^^, (2.5) 

where 5X is an infinitesimal hermitian random matrix with averages 

{6Xij) = , {6X,,5Xki) = D5t 6ik6ji (2.6) 

and that is independently distributed of S{t). The diffusion constant D will be precised 
afterwards, depending on the considered physical problem. We will only consider the unitary 
case. The brownian motion ensembles for the other cases are more involved and their precise 
definition can be found in Ref. [Q. Let p{(p, t) be the joint probability density for the angles 
(Pj of 5* at the time t, obtained after integration over the unitary matrices ^1^2? ^^1,2- This 
Brownian motion is then characterized by the Fokker-Planck equation: 
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(2.7) 



dt 4 Y V d^. 



with Fp{(p) given by Eq. (p.4|) . In the Appendix this Fokker-Planck equation is derived for 
the unitary case /3 = 2 which is our main concern. 

The details of the derivation in the Appendix show that the Fokker-Planck equation ( p.7|) 
remains valid even if the statistics of the matrix 6X is arbitrarily chosen (including the case 
of a constant matrix 6X) instead of ( |2.6| ). For this to be true, one then needs instead of 
Eq. (p.6|) the non-trivial assumption that the unitary matrices V2 and U2 that appear in the 
parametrization ( p.l|) of S{t) are for each t uniformly distributed on the space of N x N 
unitary matrices (i.e. they have always two independent N x N CUE distributions). The 
time step 6t is then determined by Eq. ( [A8| ) of the Appendix. It is interesting to note that 
this assumption concerning U2 and V2 is fulfilled at least at t = for the two applications 
concerning the crossovers 2xCUE CUE or COE CUE which are treated here or in 
Ref. respectively. We assume that this holds also for arbitrary t > and the Brownian 



motion model is applicable to rather general physical situations. 



B. Solution of the Fokker-Planck equation for arbitrary initial condition 



We want now to solve the Fokker-Planck equation ( p.7|) for an arbitrary initial condition 
p{(p,t = 0) = p{(p) and for (3 = 2. We set D = A for convenience. We proceed in a similar 
way as in Ref. and apply the so-called Sutherland transformation [p6| : 



p{ip,t) = iFpiif)r'^'piif,t) . (2.8) 
The function p{ip, t) fuUfils a Schrodinger equation with imaginary time: 

dp{<^, t) 



dt 

where 7i is a many particle Hamilton operator 

2 



np{ip,t) (2.9) 



n = -Y.(i^] +yi^) (2-10) 



with a potential V{(f) having the form 



Vi^) = - E . 2,0 ^ (si^' + {P-m- P) cos' 



J sin {2ipj 

j^k (sin' Lfj - sin' ip^j 



^ /?(/3-2) sin^(2(^,-)+sin^(2y.fc) ^ ^ 



and the constant 



1 B' 1 

Cj^ = —3'N -—N{N -1){N -2)- -(3{(3 + 2)N{N -I) . (2.12) 
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One can see that the unitary case /3 = 2 is apparently much more simpler, since the com- 
plicated many particle Hamilton operator reduces in this case to a sum of independent one 
particle Hamiltonians, i.e. Ti = J2j h{ipj) + with 



The eigenvalue problem of this operator can be solved using Legendre polynomials. The 
properly normalized eigenfunctions ipniVj) with eigenvalues En, n = 0, 1, 2, . . ., are given by 



M^j) = V (1 + sin(2y^j) P„(cos(2y;,-)) , = (1 + 2n)^ . (2.14) 

The solution of the Fokker-Planck equation can be expressed by the many particle fermionic 
Green's function G{ip,ip;t), as already done in Ref. 



via 



p{if,t)= j d''(pp{(p)G{(p,^-t) (2.15) 

where 

G(<^, t) = ^Piv)'' P{^) det{g{ip,, t)) e"^-* , (2.16) 

P('P) = n (sin^ "Pi - sin^ "Pj) 11 \/sin(2(pj) , (2.17) 

i<j j 
oo 

g{(p^, (pj; t) = Y^ iJni^Pi) M^j) e"""* ■ (2.18) 

n=0 

Eq. (|2.18|) describes the one particle Green's function, defined by 

(^^ + h{v>i)^ giVi^'Pft) =0 , g{ipi,ipj;0) = 6{ipi - <pj) . (2.19) 

For technical reasons, we will now switch to new coordinates Xj = cos{2(fj) G [—1, 1] 
which are just the arguments of the Legendre polynomials. Let p{x,t) be the probabil- 
ity density in terms of the xj taking into account the Jacobian of this transformation, i.e. 
Y\j{4:sm{2{pj))p{ip,t) — ^ p{x,t). Eqs. (|2.15| - p.l8D are then replaced by 

p{x,t) = J d^x p{x)G{x,x;t) , (2.20) 

G{x, X- 1) = ^P(5^)"' p{x) det(^?(x„ X,- 1)) e-^'-* , (2.21) 

p{x) = X{{x,-x,) , (2.22) 

i>j 
oo 

g{xi, xf,t) = Y.-{l + 2n) P^Xi) P„(x,) e"^"* . (2.23) 

n=0 ^ 



We have used in ( |2.20[ - p^23[ ) the same symbols p, p, G, g and p for objects which are related 
with the corresponding objects of (p.lSj - pTTsD via suitable transformations including the 
corresponding jacobians, and some factors Jsm{2(pj). 
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One can easily check that p{x, t) in the hmit t — > oo corresponds to the CUE distribution. 
For this, we need to consider only the first N contributions, i.e. n = 0, 1, . . . , — 1, in the 
sum of ( p.23| ). The matrix g{xi,Xj;t) is then the matrix product of two matrices P„(xj), 
Pn{xj) and a diagonal matrix with entries (| + n)e~^"*. Then the determinant of Eq. (p.21|) 
becomes proportional to the two Vandermond determinants p{x), p{x) and to the factor 
exp(— Y,n=o ^nt) which is just canceled by the contribution e"*"^* in Eq. (|2.21|) . This gives: 



lim t) oc p^(x) ■ / d^x p{x) (2.24) 

t^oo J 

which corresponds to a stationary solution equal to F2{ip), once the jacobian of the transfor- 
mation Xj = cos{2(pj) has been taken into account. The first correction to the limit ( p.24| ) is 
proportional to exp{—{eN — £N-i)t) = exp{—SNt), giving the typical time scale tc = 1/{8N) 
for the Brownian motion to reach the stationary solution. 



C. Universal parametric correlations and their limits 

The universality of the energy level parametric correlation has recently attracted a con- 
siderable interest. In mesoscopic quantum physics, this remarkable property obtained after 
an appropriate rescaling of the variables, was proven by diagrammatic methods in Ref. [§. 
An alternative derivation of this result can be directly applied to our problem: Beenakker 
and Rejaei have indeed shown how to recover this universal behavior from a Brownian 
motion model for the Hamiltonian (//-Brownian motion ensemble). This derivation can 
be easily adapted to the S'-Brownian motion ensemble and to the transmission eigenvalue 
parametric correlation. As in Refs. we do not consider here a crossover regime (i.e: 

2x CUE ^ CUE or COE ^ CUE for instance) but just a CUE ^ CUE parametric depen- 
dence, which reduces to the large t-behavior of the crossover of two decoupled N x N CUE 
towards a single 2A^ x 2A^ CUE considered in section pT[ Our discussion follows very closely 
the derivation of Ref. |Q, where a parametric dependence on a quantity X was considered, 
related to the Brownian motion time through the relation X"^ = 8Nt = t/tc- We measure t 
in units of 8A^ = in order to have a stationary distribution proportional to , and we 
use again the variables ip. [] 

Usually, one considers two types of correlation functions which are the density 
correlation function 

s{^, X, X') = y: - v^{x)) s{v' - v,{x'))) 

- E - v^ix))) y: m^' - ^.(xo)) (2.25) 

i j 

and the current correlation function 

Cicp, X, X') = Y: {(HX) MX') - V^iX)) 6icp' - MX'))) . (2.26) 



We think that they are the appropriate coordinates since they have a uniform density in the 
Hmit N ^ 00. 
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For an arbitrary linear statistic (e. g. conductance, shot noise power,...) of the form 



A{X) =Y.a{v,{X)) (2.27) 



we can calculate the correlator via 

f7r/2 i-k/2 



{SA{X)5A{X')) = r dip r dip' a{ip)a{ip')S{ip,X,ip',X') . (2.28) 
Jo Jo 

The two correlation functions ( p.25|) and ( |2.26|) are related through 



S' „ 



In subsection p.11 A| , we calculate the probability density Ri^i{x,t;y,t + r) to find the 
eigenvalue x at time t and the eigenvalue y at time t + t. For a given r and in the limit 
t oo, the particular information of the initial condition will be lost and we have the 
density correlation function for the CUE — > CUE parametric dependence: 

Sn{^,0,^',X) = lim 4sin(2(^) sin(2(^') (2.30) 
x(i?i^i(cos(2v9),t; cos(2v9'),t + r) - Ri{cos{2ip);t) Ri{cos{2ip');t + r)^ . 

The parameter X is now related with r through = r/t^ and the sinus-prefactors are 
just the Jacobians due the variable transformation x = cos(2(y9). The subscript means that 
N can have arbitrary values and that the limit N ^ oo has still not been taken. Using 
the results ( |3.25| ), ( |3.26| ) and ( |3.12| ) derived in the next section, we get directly the exact 
expression: 

SN{'P,0,ip',X) =4sin(2y?) sin(2(^') 

X ( £ + i) ^n(cos(2¥;)) P,(cos(2y.')) e-'A (2.31) 

\n=N J 

X f E + I) Pkicos{2^)) P,(cos(2(^')) e+'A 

\k=0 J 

where = (1 + 2n)^ is the eigenvalue of the one particle Hamiltonian h{ip) given by ( [^.13[ ). 
In another work [^], we have calculated the function Ri^i for the crossover COE CUE 
which differs at small t from the result ( |3.25 ) we have used here. However, in the limit 
t — *• oo both functions become identical and give the expression (|2.31|) which has in fact the 
same structure than in Eq. (5.27) of Ref. P]. One has just to replace the harmonic oscillator 
eigenfunctions used in |^ by the Legendre polynomials appearing in the S-Brownian motion. 

First, we consider the limit ^ oo and differences 6ip = ip' — ip of the order 1/A^, in 
order to use the large n expansion of the Legendre polynomials: 




P„(cos(2(^)) ~ J , cos ({2n + 1)^ - f ) (2.32) 

\' /sin(2(^) ^ ^ 
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and to replace the sums by integrals. After some smoothing over the fast oscillating terms, 
the correlation function depends only on the difference S^p, i.e. 



S[5^,X)= Wm SN{ip,Q,V + 5ip,X) (2.33) 



pI J ds J ds' cos{po'iTS 5(p) cos{poiTs' 5(p) exp (|aX^(s^ — s'^ 



where po = ~" is the average density of the variable ip and a = N/2 = Trpo/A. The expression 
( p.33| ), which is valid in the large N- and small 6(p-limit, coincides exactly with the result of 
Ref. 1^ given by a Brownian motion for the Hamiltonian H and also with a supersymmetric 
calculation for the unitary case The calculations of the current correlation function 
( p.26| ) and of the correlator p.28| ) are then straightforward and can be found in Ref. [§] . We 



note that (|2.33| ) yields an algebraic decay for the correlator of an arbitrary linear statistic 



as X ^ oo § 



(M(0)M(X))~^ , (2.34) 



which is also found pO[ in a simplified Brownian motion for the transmission eigenvalues. 
However, we cannot directly use the expression (|2.33| ) to evaluate the integrals in ( p. 28 ), 



due to the restriction 6ip <^ 1. In fact, we will see in the next section that the conductance 
correlator decays exponentially in the S-Brownian motion ensemble. In order to understand 
this behavior, we consider again the limit N oo, but now without any restriction on ip 
and (fi'. We replace in the exact expression ( |2.31| ) again the Legendre Polynomials by the 
asymptotic formula ( 2.32| ) but now we keep the discrete sums instead of replacing them 



by integrals. One can see that this step is very crucial: the variable varies between 
and 7r/2 and the correlation function should rather be expanded in a discrete Fourier-series 
instead of a continuous Fourier-integral, being responsible [|] for the algebraic decay ( ^.34p 



In the limit ^ oo at fixed ratio r/tc = X^, the main contributions comes from the 
terms with n — N <^ N and N — k <^ N so that we can linearize the energy difference 
^k — = 4(1 + k + n){k — n) ~ 8N{k — n). In addition, all fast oscillating terms of the type 
cos(4A^(- ■■) + ■■■) are omited. Then, one of the sums is easily done and we obtain for the 
density correlation function the expression 

lim Sn{^, 0, X) = S{ip + if', X) + S{ip - ip', X) (2.35) 

N^oo 

where the function S{ip,X) is given as a Fourier- series 

oo oo 

Siifi, X) = 2J2 Sn{X) cos(2n¥?) = ^ Sn{X) e^'""^ (2.36) 

71=1 n=— OO 

with Fourier-coefficients 



5.(X) = f^e-H^' , (3 = 1,2, A . (2.37) 



Our derivation accounts of course only for the unitary case (3 = 2 but Eq. ( p.37|) is valid for 



the orthogonal (/? = 1) and symplectic case (/? = 4) too [^. The derivation for all (3 can be 
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done by an asymptotic expansion of the Fokker-Planck equation ||3J] similarly as in section 
4 of Ref. P and Eq. ( p.37|) is indeed the discrete version of Eq. (4.14) in Ref. [^. 



We are now able to calculate the correlator for the linear statistics ( p.27| ). If we expand 
the function a{ip) in a Fourier-series 



oo 



we find 



a{^)= ' an = a_n , (2.38) 



OO O OO 

(M(0)M(X))=7r2^a^5„(X) = -^a>e--^^ . (2.39) 

n=l P n=l 



In ( 2.38| ), we have extended a{(p) to an even function on the interval — 7r/2 < < vr/2 



so that the Fourier-coefficients are even in n (One should imagine that a is a well defined 
function of x = cos(2(y9).). 

As an example we consider the conductance g = A with a{ip) = siia^{(p), i.e. Oq = 1/2, 
a-ti = —1/4 and a„ = if \n\ > 2, 

69{X)) = ^ e-^^ = ^ e-^/*= (2.40) 

which is (in the limit t —>■ oo) consistent with Eq. (|3.30|) of the next section. 

The current correlation function is calculated from ( p.29| ) and the assumption that S 
and C depend only on the difference X' — X. We obtain 

C{^, 0, X) = -C{^ + X) + - X) (2.41) 

where C{ip,X) has the Fourier-coefficients 

a(X) = -^^5„(X) = ^(l-2HX^)e-H^^ , n^O . (2.42) 

The coefficient Co{X) is not determined and not needed due to ( |2.41|) . The modified sign in 
the first contribution in ( p.41| ) is important because the current correlation function is odd 
as a function of ip (or if') if ip' (or ip) is fixed. 

In summary, we have found that the parametric correlations of the S'-Brownian motion 
are only locally (on a range 6(p = (p' — (p <^ 1) identical with the universal correlations 
found by a i7-Brownian motion |Q] or by a microscopic approach in the framework of the 
zero-dimensional super symmetric non linear cr- mo del 0. This agreement does not hold if 
we consider the whole range of ip, ip'. In particular, quantities integrated over the whole 
spectrum like the correlator of a linear statistic must be treated with more care and indeed 
the conductance correlator shows an exponential instead of an algebraic decay, in the S- 
Brownian motion ensemble. 
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III. APPLICATION TO TWO WEAKLY COUPLED CHAOTIC DOTS. 



The S'-Brownian motion ensemble with CUE stationary hmit can be apphed to a system 
made by two quantum dots in series under an apphed magnetic field, for describing how its 
transmission eigenvalue distribution varies as a function of the strength of their coupling. 
We assume that the two cavities are coupled between them by a wire with n conducting 
channels (see Fig. 1), and each of them to one electron reservoir through n) channels. 
The underlying classical dynamics is supposed to be fully chaotic in each dots. 

For this particular application, we need to express the fictitious time t of the Brownian 
motion, which has been introduced in a rather abstract mathematical way, in terms of the 
physical parameters and n. In Ref. it was shown that a closed system where the 



classical motion remains (nearly) in several separated chaotic regions of the classical phase 
space can be modelled by a random matrix approach similar to Dyson's Gaussian Brownian 
motion ensemble discussed in Ref. @]. The time parameter of the latter can then be related 
with the particle fiux between the nearly separated regions of phase space. In principle, one 
can relate pBf the if-Brownian motion (for a closed system) with the S'-Brownian motion (for 
an open system) by some microscopic assumptions about the S'(if)-dependence and then 



apply the results of Ref. [§2| in order to obtain the physical meaning of the S'-Brownian 
motion time. 

For the scope of this paper, we will only give a simple heuristic derivation. Therefore, 
we assume that the (A^ + n) -dimensional S'-matrix of each cavity has a CUE distribution 
when an external parameter is varied (Fermi energy, magnetic field, shape of the dot...) and 
has the form shown in Eq. ( |1 . 1|) , excepted that r is x iV-matrix whereas the coupling 
gives rise to a n x A^-matrix t. When one electron is injected from a reservoir into one of the 
cavities by the z-th of the A^ incoming channels of the corresponding lead, one can assume 
that the probability that it goes directly back to reservoir is just 

N 

and that the probability to go through the n channel wire into the other cavity is 

Ptrans = ^ / = nr , = ^ " Pout ■ (3.2) 

j=i i\ +n 

This results from the assumption that the chaotic motion inside the dot uniformly explores 
the dot boundaries. Furthermore, one can assume that the typical time for the electron to 
stay in one cavity is very long compared to the time needed to cross the connection, i.e. the 
electron looses completely the memory from which reservoir it was injected to the actual 
cavity. In this case, the total probability ptotai for the electron to be transmitted from one 
reservoir to the other is given by ptotai = P trans (p out +P trans P total), a relation which accounts 
for all type of processes where the electron can move between the cavities an arbitrary 
number of times before it leaves the system into the other reservoir. The solution of this 
equation is given by 

Ptrans ^ fo o\ 
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Then, the average conductance can be estimated by {g) = A^p total which leads with the 
exact result Eq. ( p. 27] ) of Section III A to the identification 



t:^J-ln(l+'^) (3.4) 

where t the Brownian motion time. The prefactor in (|3.4|) accounts for the typical time scale 
tc = 1/ {8N) which is found in the end of Section p.1 lj| . In the weak coupling limit n -C A^, 
we have the simple proportionality t/tc — {2n)/N. 



A. Crossover of two CUEs to one CUE for the transmission eigenvalues 



In the parametrization (|2.1| ), the case of two independent CUEs of dimension A^ cor- 
responds to the requirement that all Tj = or Xj = 1 respectively. We have therefore 
p{x) = — Xj) for the initial condition. The evaluation of (|2.20| ) leads to the limit 

p{x, t) = lim^^^i G{x, x; t). We note that the Legendre polynomials are particular Hyperge- 
ometric functions, i.e. P„(x) = F{—n, + 1; 1; i(l — x)) which can be expanded as 

Pnix) = j:qkUn+'^r) (x-l)' (3.5) 



with the polynomials 

^^M = ^n(«-(-+i)^) (3-6) 

i/=0 



of degree k. We now expand the one particle propagator g{xi, Xj] t) with respect to (1 — Xj) 
up to terms (1 — Xj)^~^. Again we can identify a matrix product in the determinant of Eq. 
( p.20| ) which yields a factor p{x) just cancelling the inverse contribution in ( p.20|) . For the 



further calculation we use the notation Pn{x,t) = Pn{x) e*^"* and find from ( |2.2(J| ), ( |2.21| ) 



p{x, t) oc det(P„(xi, t)) det{gk{xi, t)) (3.7) 

with 

oo 

gkix,t) = Y,{n + l)qk {{n + i)^) P„(x) e"^"* . (3.8) 

n=0 

Because of the determinant in Eq. (|3.7|) , we note that one can replace in Eq. ( p.SD the 
polynomials qk{u), k = 0,1, . . . , N — 1 by an arbitrary set of linear independent polynomials 
of (maximal) degree A^ — 1. Such a substitution affects only the normalization constant. 
In the following, we will use two different substitutions for two different purposes: first 
to calculate exactly all kind of correlation functions and second to determine the effective 
interaction between the xj (or (pj) variables. We begin with the calculation of the correlation 
functions defined by 

A^! f 

Rk{xi,..., Xk, t)= _j^y j dxk+i ■ ■ ■ dxN p{x, t) , k=l,2,...,N . (3.9) 
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We put ui = {I + = jEi, I = 0,1, . . . , N — 1 and replace the polynomials qk{u) by the 
Lagrangian interpolation polynomials Lk{u) for the Uk, i.e. 



Lk(u) 



N-l 

n 



Uk - Ui 



Lkiui) = 6, 



kl 



(3.10) 



where I, k = 0,1, . . . , N — 1. We find now for the joint probability density the result 

1 



with 



p{x,t) = -^det{Pn{xi,t)) det{hk{xi,t)) 



hk{x, t) = {k + i) Pkix) e^^^* + E + I) Lk {{n + ^f) P„(a;) e"^"* 

n=N 



(3.11) 



(3.12) 



The normalization in ( |3.11| ) follows directly from the obvious property 

j ^dxPn{x,t)hk{x,t) = 5nk ■ 

Actually, this "orthogonality" relation enables us to calculate all correlation functions ( p^ 
We define the function 



(3.13) 



N-l 



Km{x, y;t)= E Pn{x, t) hn{y, t) 



(3.14) 



n=0 



and find the following three properties 



p{.x, ^) = ^ det {KN^Xi, Xj] t) i<Lj<N) 



dx Kn{x, x;t) = N 



J ^ dy Kn{x, y; t) K^iy, z; t) = Kn{x, z; t) 



(3.15) 
(3.16) 
(3.17) 



From these properties and a well known theorem in the theory of random matrices PJ37 
(cp. Theorem 5.2.1. of Ref. [Q), we find directly the result for the correlation functions 



Rk{xi, ...,Xk;t) = det {KN{xi, xf, t) i<i,j<k) ■ 

For the cases k = 1,2, we have 

Ri{x; t) = Kn{x, x; t) , 
R2{x, y; t) = Ri{x; t) Ri{y; t) - Kn{x, y; t) KN{y, x; t) . 



(3.18) 



(3.19) 
(3.20) 



For practical applications, we give the more explicit expression for (|3.14|) 



N-l 



N-l oo 



K^{x,y-t)= {l + 2k)Pk{x)Pk{y)^Y. Y.i^ + \) DknPk{x) Pn{y) e^"--'-^' (3.21) 

fc=0 fc=0 n=N 
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with coefficients 



I?.. - L, ((n + ,) j - Jl^^ + ; + (3-22) 

(_l)JV-i-fe(^^ jY)! (2A; + 1)(A^- /t) 



(ri-iV)!(A; + iV)!(Ar-fc)! (n + + l)(n - A;) 
One might also consider correlation functions at two different times, i.e. 

Rk,i{xi, . . . ,Xk,t;yi, . . . ,yi]t + r) 



(3.23) 



(A^- A;)!(A^-/)! 

X J dxk+i - ■ -dxN J dyi+i - ■ -dyN pix,t) G{x,y;T) . (3.24) 

Without going into technical details, we mention that one can calculate the function 
Ri^i{x,t; y,t) with the technic of functional derivatives [Q. The result is 

RiA{x,t;y,t + t) = Ri{x;t) Ri{y;t + T) + (^g{x,y;T) - KN{x,t;y;t + T)'j KN{y,t + t; x,t) 

(3.25) 

where 

Kn{x, ti; y\ ti) = J2 ti) hkiy, (3.26) 

fc=0 

is a generalization of ( p. 141) for different times. The limit r ^ is given by i?i^i(x, t; ?/, t) = 
R2{x,y;t) + 6{x - y) Ri{x). 

The knowledge of the correlations functions ( p.l9| , |3.25| ) enables us to calculate exact 
expressions for the average of the conductance g = J2jTj = |Sj(l ~ a^j) and the auto 
correlation of its variance Sg = g — (g). Using the orthogonality of the Legendre polynomials 
and their recursion relation we find 

(^7(t)) = f (1 - e-«^*) , (3.27) 

mt) 6g{t + r)) = ^ e-- |— ^ + i±\ e^^-^^)* 



4 UA^2 _i 2N+1 

+Z_^e-«(2^-i)* -1-e-i^^*} . (3.28) 

Of particular interest is the large limit. In order to obtain a non trivial crossover, we 
have to measure the time variable in units of the typical crossover scale tc = 1/ {8N) and to 
keep the ratio t/tc fixed in the limit ^ 1. We find then up to corrections of order 1/A^ 

(^W) = f (1 - e-*/*^) , (3.29) 

m)5git + r)) = ^e-^/^^ ll+(^-f--f-l\e--^/A . (3.30) 
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B. Effective pairwise interaction in the crossover regime 

We want evaluate ( |3.7D in a suitable way for having the effective pairwise interaction 
between the Xj- or y^j-variables. We now replace in Eq. ( |3.8|) the polynomials qk (^{n + 
by 4^(n + D^*"' = and find 



with 



p{x,t) oc det{Pn{xi,t)) det{gk{xi,t)) (3.31) 



oo / \ 

hix, t) = X: (n + i) et Pn{x) e-^"* = (-1)'= ^-j go{x, t) 



f7o(x,t(l-s))^^^ . (3.32) 



The quantity go{x,t) can be related to the heat kernel on the unit sphere 5*2 with initial 
condition at the origin, which yields a useful approximation in the limit t <^ 1 



From the expansion 



^7o(cos(2^),t)^le-^'/(^*) . (3.33) 



1 _^^/n_o^ x-^ Lhix 



■ e 



xs/{l-s) _ sr^ ^k\X) 



E (3.34) 



l-s k\ 
(here L^^x) denote the Laguerre polynomials) and ( |3.33| ), ( |3.32| ) we find 

^.(cos(2^),t) = ^e-'^V(^*)L,(^^^ , t«l . (3.35) 

This results leads to a probability density (now for the variables (pj ) 

p(^, t) ^ C{t) n (sin(2v.,) e-^?/(^*)) 1] { (sm^ - sin^ (^^ - ^^2) } (3.36) 

where C{t) is a normalization constant which depends on t. If we write ( |3.36| ) under the 
form of a "Gibbs factor" p{(p,t) ~ exp[—P{J2i<j u{(pi, (pj) + J2iV{'^i,t))], we find that the 
fictitious corresponding Hamiltonian has the pairwise interaction and the one body confining 
potential: 

u{Lpi, Lfj) = (in I sin^ - sin^ ipj \ + In - V^jl) , (3.37) 
r(¥.„t) = |i-^ln|sin(2^,)| . (3.38) 

This can be simplified when t is very small, i.e. t <^tc = 1/{SN) since sin^ ipi ~ ip] in this 
limit. One gets a joint probability distribution (|3.36|) identical to the so-called "Laguerre 
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ensemble" of Ref. [|38|, with a logarithmic repulsion between the ip] (see also Ref. ||3^). For 
N ^ 1, the density reduces to a quarter circle law: 

1 



Ri{^;t)^—^mN-ip^ (3.39) 

if the argument of the square root is not negative (otherwise Ri{(p; t) = 0). Using Lp^ c::iT = 
cosh~^ z/, one gets 

t) = ^^^i?i(cosh-^ z/; t) (3.40) 
cosh V 

for the variable v G [0, cxo]. 

One can note that the exact pairwise interaction for the (pj variables ( p.37|) differs for 



from the standard RMT interaction given by ( p.4|) for (3 = 2. However, this 
concerns essentially the quantitative difference between sin^ ip and ip"^ whereas for the quasi 
one-dimensional disordered wire considered in Ref. this effect was much more stronger 
due to the appearance of sinh^(. . .) instead of sin^(. . .). 

In the sections 0, we will return to this point and compare the results Eqs. ( p.l9D , 



( ^.4(J| ) with the corresponding expressions obtained for disordered lines, assuming the M- 
Brownian motion ensemble. 



IV. M BROWNIAN MOTION ENSEMBLE FOR THE TRANSMISSION 

EIGENVALUES 

The S-Brownian motion ensemble was built on the evolution law S{t + 6t) = 
S(t) exp{i6X) where 6X was an infinitesimal hermitian random matrix (Eq. ( |2.5|) ). A 
Fokker-Planck equation was deduced, giving the t-dependence of the transmission eigen- 
values. A different Fokker-Planck equation is more natural when one builds a quasi-one 
dimensional wire, adding in series small diffusive blocs p7|-p0|. The evolution law of the 



matrices u and v ( "angular part" ) and of the matrix T ( "radial part" ) defined in Eq. (ETT 



can indeed be deduced from the multiplicative transfer matrix M, expressed [E3 in the same 



parametrization. The law of M describing the addition in series of the scatterers is now: 

M{t + 6t) = M{t)M{6t) . (4.1) 

The Brownian motion time t is in this case related to the number of scatterers in series 
(i. e.: the length L of the wire). M can be parametrized with the same coordinates used 
for S in Eq. 2A, as detailed in Ref. |jl5|] for instance. Assuming that the scatterers are 



isotropically distributed, the length dependence of the transmission eigenvalues for the M- 
Brownian motion ensemble is given by a Fokker-planck equation [p^^-|30[], which presents 



interesting analogies and differences with Eq. (|2.7| ). The comparison of those two Brownian 
motions, assuming two different evolution laws, will enable us to describe in details the 
similarities and the differences between quantum transmission through two weakly coupled 
cavities (S'-ensemble) and through a disordered wire (M-ensemble). Before this, having 
used the Beenakker-Rejaei method for solving the S'-Brownian motion ensemble, we have 
been able to improve this method and we show in return how to complete the solution of the 
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transmission problem through a disordered wire. Indeed, we explain how to calculate exactly 
the length dependence of n order correlation functions of the transmission eigenvalues, for 
arbitrary values of n. This allowes us to establish an important equivalence between the M- 
Brownian motion ensemble and a more microscopic approach using the non linear cr-model. 



A. Exact calculation of the correlation functions for (3 = 2 

We start with the exact result given in Ref. [^] for the joint probability distribution of 
the variables \j = (1 — Tj)/Tj = sinh xj which satisfies the Fokker-Planck equation of the 
M-Brownian motion ensemble: 

p(A, t) oc lliXi - Xj) det{gm{Xj, t)) e^^* , (4.2) 

i>j 

with 

gm{Kt) = r dk{erltanh{f) P.,^,^,^{l + 2X)e-''' , (4.3) 

^ 

where Cat is just given by (|2.12| ) for (3 = 2. We underline that the changed sign in the 
exponential factor of (|4.2| ) if compared with (|2.21|) is correct. The variable t is now t = 
L/{4Nl) = L/{2C,) {I is the elastic mean free path, L is the length of the conductor, and 
^ = 2NI is the quasi one-dimensional localization length for the unitary case). The index m 
in ( [4.3| ) takes the values m = 0, 1, . . . , — 1. For later use we mention that the generalized 
Legendre function can be expressed as a Hypergeometric function via 

Pi(.fc_i)(l + 2A) = F(i-z|,i + z|;l;-A) . (4.4) 

This function is an eigenfunctions of the differential operator 

with eigenvalue fc^. The operator D is of course related with the one particle Hamiltonian 
Ho used in Ref. by a suitable transformation. We can now state the following properties 
of gm{X,t) 

gmiKt) = D^goiX,t) , (4.6) 
^^go{X,t) = -Dgo{X,t) , (4.7) 
goiX,0) = 6iX) . (4.8) 

Eqs. ( [4. 61 ), ( [4.7] ) follow directly from ( (4.31 ). Eq. ( [4. 81 ) reflects the initial condition for the 
one particle Green's function and was already given in |^, where the expression ( [4.2|) 



was found for the one channel case. 

We are now interested in integrals of the form 



it) = / dX P„(l + 2A)(7n^(A,^) (4.9) 

Jo 
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where f„(- ■ ■) again denotes the Legendre polynomial of degree n. We consider first the 
case m = 0. From (|4.8|) we find cirio(O) = -Pn,(l) = 1 and from ( |4.7D we get the differential 
equation 

d r°° 

Tranoit) = - dX P„(l + 2A) D go{X, t) 
at Jo 

POO 

= - dX go{X, t) D P„,(l + 2A) = (1 + 2nf ano{t) = Sn a„o(t) (4.10) 
Jo 

which has the solution a„o(^) = exp(e„i(:). In Eq. (^4.10|) we have used that P„(l + 2 A) is an 
eigenfunction of D with eigenvalue — = — (1 + 2nY , a property following directly from 
Legendre's differential equation. From Eqs. (|4.6| ), ( [4.9|) we get by a similar calculation the 
general result 

anM = {-euT (^no{t) = {-SnT e'-' . (4.11) 

As in Sec. [Ill A| we can replace the exponential factors (fc^)™ in Eq. ( |4.3| ) by an arbitrary 
set of linear independent Polynomials of (maximal) degree — 1 in the variable k'^ because 
the determinant in Eq. ( [4. 21 ) yields only a constant factor after such a transformation. We 
now choose the replacement 

{k^r ^ L„,(-\k') (4.12) 



where Lm{u) are just the Lagrangian interpolation polynomials already given by Eq. ( p. 101 ). 
We introduce the following notations 

Q^{X,t) = Pn{l + 2X)e-'"' , (4.13) 

and 

hr,,{X,t) = dkLrn{-\k'') | tauh (f ) Pi (1 + 2A)e-^'* . (4.14) 
From Eq. ( [4.11| ) we get directly for n, m = 0, 1, . . . , — 1 

dX Qn{X,t)hm{X,t) = Lra{\en) = Snm ■ (4.15) 

The joint probability distribution ([4.2[) can then be written in the form 



p{X,t) = ^ det{Qn{Xj,t)) det{h^{Xj,t)) = ^ det (A'Ar(Ai, A^; t) i<ij<Ar) (4.16) 

with 

N-l 

K^iX, X;t)=Y. Qm{X, t) hm{X, t) . (4.17) 

m=0 



Similarly as in Sec. |III A| we find from ( [4.15|) 



POO 

/ dX KN{X,X]t) = N (4.18) 
Jo 
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and 

/ dfi KNiX,i^;t)KNii^,X;t) = KNiX,X;t) . (4.19) 
These properties lead as in Eq. ( |3.18|) to the m-point correlation functions 

Rm{Xi, ...,Xra;t)= dct (/C7v(Ai, Xj] t) l<i,j<m) (4.20) 

where Kn{X, A; t) can also be expressed by the more explicit formula 

KNiX,X;t) = /"rfA;|tanh(f )Pi(,^_i)(l + 2A)e-'='* 
Jo ^ 

Af-1 

^L^(-ip)p„(l + 2A)e-^-* . (4.21) 



Af-l 

X 

m=0 



B. Equivalence with the microscopic non Unear cr-model approach 

As a first application we calculate an exact expression for the average conductance, when 
the length of the wire increases. This expression is valid from the disordered conductor 
towards the Anderson insulator. Since the conductance is a linear statistic of the Xi, g = 
X)i(l + Aj)"^, we have to the evaluate the corresponding integral over the density: 

{g)= dX-—-J2QmiX,t)hUX,t) . (4.22) 

Since Qm{X, t) is a polynomial of degree m in the variable 1 + A, we can write 

1 1 f — 11™ 

Qm(A,t) = — — g^(-l,t) + w(A) =^—-l-e-'-' + rrn-i{X) (4.23) 



l + A"'^ ' ' l + A"'^ ' ' ^ ' 1 + A 

where r^-iiX) is a polynomial of degree m — 1 which does not contribute in the integral 
( [4.22| ) due to the quasi orthogonality relation ( |4.15| ). We find therefore 



{g) = ^ e-(^-+'=')* (-1)- Lm (-|A;2) I tanh (f ) ■ I{k) (4.24) 

m=0 

where I{k) stands for the integral 

m = l ^F{l-^|,l + ^|■A■,-X) . (4.25) 



We use the transformation formula 15.3.4 of Ref. [^] for the Hypergeometric function in 
( [4.25| ) and expand afterwards the latter in its power series. The integral over A can then be 
done and gives the sum 



m = E 1 L 'In (4-26) 



„=0 2+^-^f V ^ 
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where denotes the generahzed Binomial coefficient which is defined for arbitrary complex 
values of a. The evaluation of the sum ( ^.26| ) yields 



du {I —U) 2+*2 -U 2 *2 



r(i + .|) r(i-.|) 



(4.27) 



We now use the expression (|3.10| ) for the polynominals 
conductance the expression 



and find for the averaged 



roo 

{g) = 2j2 dk ^-{{i+2m?+k')t k tanh ( 

■m=0 



2m + 1 



P + (1 + 2m)2 



a{N, m, k) 



(4.28) 



where we have introduced the coefficient 
A;2 + (1 + 2m)2 



a{N, m, k) 



4(2m+ 1) 

r(iv + | + .|) t{n + \- 

V{N -m) r(A^ + m + 1) 



(4.29) 
(4.30) 



The result ( [4.28| ) is exact for all values of and t = L/{2^). In the limit N ^ oo the coef- 
ficient a{N, m, k) is just 1 and the expression ([4.28|) becomes identical with the microscopic 
result of Zirnbauer et al |^5| based on the supersymmetric non linear sigma model, in the 
unitary case. 

A similar calculation yields the second moment of g 

1 /-oo 

(^2) = -Y. dk e-((i+2-)'+'=')* k tanh (f ) (2m + 1) a{N, m, k) . (4.31) 

We omit the details, since this result can also be derived in a more direct way by the general 
identity 



4 dt 



(4.32) 



which follows directly from the original Fokker-Planck equation |4^ for the case (3 = 2. 

This gives the proof that the microscopic supersymmetric approach and the M-Brownian 
motion ensemble of Ref. [0,^] are equivalent in the limit N —>■ oo, for arbitrary values of 
L, at least as far as the average and the variance of the conductance are concerned. 



C. Quasi- Id localized limit. 



In the localized limit t ^ 1 the integral in Eq. ( [4.14|) can be evaluated by a saddle 
point approximation. In order to do this, we have to consider the oscillatory behavior of the 
Legendre function for large values of the variable A. We use now the transformation formula 
15.3.8 of H and express the /c-integration measure as a product of Gamma functions. 
Instead of Eq. ( [4. 14|) we get then the modified expression 
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hm{\ t) = ^f^ dk e-^'* (1 + A)-^^t H {k, j^) (4.33) 
with the abbreviation 

H {k, ^) = F{^ + ^|,l + ^\■ 1 + ^k■ ^) (-iP) . (4.34) 

It is more convenient to replace A by the new variable z = ^ ln(l + A). A simple substitution 
then leads to the expression 

ft„.(e^= - M) ^ e-^V..,-. _1_ £ „ ^ _ - . j ,4.35) 

which is still exact. We mention that this expression is rather useful for an efficient numerical 
evaluation since the difficult oscillatory behavior of the Legendre functions has been taken 
into account by a shift of the integration path into the complex plane. In the limits t 3> 1, 
2 3> 1, the integral can be done and we obtain up to corrections of the order 1/t or e~^^ 
respectively 

hm{e'' - l,t) - ^e-^^/(^*)-^ ^^""^IXP' • (4-36) 

^ ' r(§) V4V2t^ J ^ ' 



The corresponding limit for the Qm{^,t) in Eq. ( [4.13| ) is 



Qm{e^' - 1, t) ^ e-(i+2'")'* e^^"^ . (4.37) 

— I Z 



The density Ri{z\ t) and the two point function R2{z, z] t) for the variable z can be obtained 
from Eq. ([4.20|) and Eq. ( [4.17|) . The result (containing additional factors 2e^^ or 2e^^ due 



to the jacobian of the variable transformation A — >^ 2;) is 

N-1 _ N-1 

R^{z-t)r:^Y.Gm{z-t) , i?2(^,5;t)~ ^ Gm{z-t)G^{z-t) (4.38) 

r)i=0 m,m,=0,Tny^m, 

where we have used the abbreviation 

GUz; t) = exp (-iiz- 2t{l + 2m)f) (4.39) 



for the Gaussian distribution. We mention that for (|4.38|) we have replaced the argument ^ 
of the Gamma functions and of the polynomial Lm(- ■ ■) in Eq. ( [4.36|) by the most probable 
values ^ = 1 + 2m or 1 + 2m due to the Gaussian factors. We recover the simple sum of 
Gaussian distributions with mean values 2t{\ + 2m) and a variance 2t already obtained by 



a direct simplification of the Fokker Planck equation valid in the localized limit pT| , |30| , |31| . 
The above calculation shows that the general expressions ( ^4.2U| ), ( [4.21j ) of Sec. IV A are 
indeed consistent with known results concerning the quasi-lc? insulators. 
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V. WEAKLY COUPLED CHAOTIC CAVITIES VERSUS DISORDERED LINES 



We describe now the differences between transmission through two balhstic cavities cou- 
pled by a narrow constriction (section |T|) and transmission through an homogeneously 
disordered wire of constant transverse section (section [TVl ) . This can be of practical interest 
for many purposes which are not restricted only to mesoscopic electrical conductances, but 
which could be relevant for other systems as wave guide communication lines. One can 
imagine that the signal from an antenna to a receiver is anomalously weak and that the 
radio-engineering problem is to know if the weakness of the signal is due to an homoge- 
neous deterioration of the transmission line or only to a local deterioration resulting from 
an accidental narrow constriction. 

Characteristic behaviors of the S-Brownian motion ensemble are first given, characteriz- 
ing the coupled ballistic dots under a small applied magnetic field. In Fig. 2, the density of 
the variable v = (arsinh(A))^/^ is shown for the one channel case {N = 1) at four different 
values of t/tc = 10.0, 1.0 ,0.1 ,0.01 (tc = 1/{8N) is the typical time scale introduced in 
Section |11B|) . The variables A = (1 — T)/T = {l + x)/{l—x) are the radial parameters of the 
transfer matrix, which are usual in random transfer matrix theory |r5| , p8| , p| . The curves 



are calculated from Eqs. (|3.19|) , (|3.21| ). The figure 3 gives this density for the same values of 



t but for N = 5 channels. In addition for t/tc = 0.1, 0.01 the asymptotic expression ( ^.40 ) 



that arises from the Laguerre approximation is shown, too. The latter is valid if ^ 1 and 
t -C tc- In both figures one can see that the typical values for u increase (on a logarithmic 
scale) if t decrease. One can identify a rather sharp minimum value z/min for z/ if t <^ A 
simple estimation from Eq. ( p.39| ) yields 

z^min =^ -^In f-^) • (5-1) 



2 V2t„ 

For values u > u^i^ the distribution is rather broad and the overall form is (nearly) inde- 
pendent of (apart from the case t ^ tc). A reduction of the time leads essentially to 
a translation of the density. For finite values of N (see = 5), one can see typical but 
rather small oscillations around the limiting distribution for — > oo. The value t/tc = 10.0 
corresponds within the accuracy of the plots to the stationary limit t/tc oo, i.e. the sin- 
gle CUE behavior. For this case the position of the first maximum decreases and becomes 
sharper when A^ increases. This gives rise to a conductance g ~ A^/2 (cp. Eq. ( p. 27] )). 



The corresponding behaviors of the M-Brownian motion ensemble exhibit essential dif- 
ferences from the previous case. Fig. 4 shows the i/-density for the disordered wire with 
A^ = 5 channels in three different regimes, i.e. the localized regime {L/^ = 10.0), the 
crossover regime {L/^ = 1.0) and the metallic regime {L/C, = 0.2). Figs. 5a-c contain the 
same density curves in a more readable scale than in Fig. 4. The gaussian approximation 
( ^M ) (dotted line) which fits very well the exact density at L/C, = 10.0 is indicated in Fig. 
5c. 0. The results shown in the Figs. 4, 5 were obtained by a numerical evaluation of the 



^The variable z used in Section EV C is only approximately equal to ( i.e. z = ln(cosh z^) ~ 
when v ^ 1) 
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integral (|4.35| ) except for the limits t <^ 1, u < 1 where the integral ([4.14|) is better suited. 



The substitution z = ln(cosh u) was of course exactly considered. 

One can see that the //-density of the disordered line differs essentially from the density 
shown in Fig. 3 for the two weakly coupled cavities. For small transmission ( i.e. the localized 
limit), the i/-density is mainly concentrated around maximal values i'max,m. = (1 + 2m)L/^ 
(m = 0, 1, . . . , A^ — 1) with a relative variance ((z/ — z/max,™)^)^^^/ t'max.m ~ vi/L- Between the 



maxima, the density nearly vanishes. This effect can be seen in Fig. 5c and is still increased 
if the ratio L/^ is increased. But in the crossover regime L ~ ^ (Fig- 5b), the minima are 
no longer negligible compared to the maxima. In Fig. 5a ( metallic regime), the amplitude 
of the finite oscillations are small but still clearly observable. Further calculations for 
different values of show that the amplitude of the oscillations remains unchanged if N 
is increased for a fixed ratio L/S,- In order to reach a metallic regime where the density is 
really a smooth function of one has to increase N to very high values and at the same 
time to decrease L/^ (within the limit L/^ 3> 1/-^)- 

Finally, we compare {g) and of {5g'^) between the two cases in the limit of small trans- 
mission, i.e. t ^t(. for the 2xCUE CUE crossover or L ^ ^ for the disordered line. For 
the two weakly coupled cavities, we get in the lowest order oit/tc from Eqs. ( |3.27| ), ( |3.28| ) 

t 

{g{t))^--^n , (5.2) 

where n ~ {N/2) ■ [t/tc] is the effective number of coupling channels introduced in Section 
|. These results lead to the limit 

M| = i_ (5.4) 

for the relative variance of g. The corresponding expressions for the quasi one-dimensional 
line in the localized limit L ^ ^ can be obtained from Eqs. ( |4.28|) , ( [4.32|) by an evaluation 
of the fc-integral in a saddle point approximation 

,3/= r(iv + i)^/Ly''= 



. (5.6) 
The relative variance of g is now infinite and the infinite length limit is characterized by 

L^oc {g{L)} 4 

The comparison with ( p.4| ) shows again a qualitative difference in the statistics of the con- 
ductance between the two cases. Eqs. ( ^.5|) and ( ^.6|) illustrate the well known fact that the 
average of the conductance does not significantly characterize the conductance properties in 
the localized limit. On the other hand, the limit ( |5.4| ) for the two weakly coupled ballistic 
cavities is a small quantity for a large number of channels. 
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VI. SUMMARY 



The S and M Brownian motion ensembles which we have used are characterized by 
two Fokker-Planck equations expressing isotropic diffusion equations on the compact space 
of unitary matrices (5*) or on a non compact symmetric space of pseudo-unitary transfer 

ii 



matrices 



respectively. The differential operators which appear in the Fokker-Planck 
equations are just the radial part of the Laplace-Beltrami operators of the corresponding 
spaces, describing free isotropic diffusion on those curved spaces. In the case of the M- 
Brownian motion ensemble, the variable u is the good coordinate since it measures in some 
way the geodesic distance on the space of pseudo-unitary matrices (cp. Ref. Fig. 4 



shows very clearly the free diffusion of z/. Since the range of u is non compact the rather 
strong level repulsion gives well separated maxima. For the compact space of unitary ma- 
trices, it is the variable if used in the Sections P A| and [II B| which measures the geodesic 
distance. The range for ip is compact and diffusion leads after a long time to an essentially 
uniform distribution. The level repulsion only causes some small oscillations around this 
distribution. 

Apart from the (important) differences which we have reviewed, we can also notice nice 
mathematical similarities between the two Fokker-Planck approaches. Many formulas and 
results of Section y coincide with the results of Ref. and Section |V], after replacing 
the variables ipj by the uj (= Xj in [Q) and sin, cos by sinh, cosh. This is partly related 
to the similarity of the parametrization with the corresponding parametrization of the 
transfer matrix (with the substitution T —>■ —A). This occurs only for the unitary case 
{(3 = 2). For the orthogonal and unitary cases (/5 = 1,4) this similarity does not exist, since 
the unitary matrices ^1^25 ^1,2 in ( |2.1| ) are related in a different way than for the transfer 



matrix ||23| by time reversal symmetry. 

Our results (Eqs. ( p.2| )-( p^ are valid for an arbitrary number A^ = l,2,3,...of channels. 
In the case of a single channel, the quotient of the Gamma functions takes the value 7r/4 and 
Eq. ( |5.5| ) becomes identical with the result of Ref. |Q obtained using a supervector model 
for the one- dimensional white noise potential at small disorder. Indeed, it was known that 
the one-dimensional white noise model for small disorder leads to the same Fokker- 

Planck equation as the one of Refs. p7t^p8[ for one channel. In the opposite limit N ^ 00, 
the Gamma function quotient reduces to one and we recover the corresponding expression 
of Ref. 1^ obtained by the supersymmetric non linear a- model: the expressions ( [4. 281) , 
( [4.31| ) for (g) and {g"^) coinciding exactly for all length scales < L/C, < 00 with those of 
Ref. [^. The non linear sigma model for quasi one- dimensional disordered conductors can 
at least be derived assuming three different microscopic Hamiltonians, i.e. the white noise 
potential Wegner's orbital model [Q, and a random banded matrix Hamiltonian ^ 



It is then very remarkable that in the unitary case, these microscopic approaches become 
equivalent to the M-Brownian motion ensemble of Ref. PB|. This behavior is apparently not 



true for the symplectic case, since the non-linear sigma model predicts finite values of 
(g) and {Sg"^) in the localized regime, while the Fokker-Planck equation gives an exponential 
decrease of the typical conductance for each /? = 1,2,4 ||30|. In principle, it could be still 
possible that in the large N limit the average limL^oo{g{L)) takes a non vanishing value 
because of some very subtle effects concerning the tails of the normal distribution of lii{g). 
But even this cannot happen, as it can be easily seen from the generalization [H0| of Eq. 
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( [4.32| ) for arbitrary P 



d{9l 
dL 




{92))] . (6.1) 



Here ~ denotes the localization length and g2 is just the trace of the squared trans- 
mission matrix (tt^)^, i.e. g2 = + From the inequality < g2 < one finds 



directly that in the limit L 00 the behavior (g) = const. 7^ (as predicted in Ref. ||35 
for /3 = 4) is not possible since {g"^) must vanish, which yields directly (g) = 0. This contra- 
diction when (3 = 4 between the two approaches is surely worth to be considered in future 
work. 

Though the two Fokker-Planck approaches are very similar from a mathematical point 
of view, they nevertheless describe two different physical situations where the Fokker-Planck 
time controls the crossover between systems with a high and a low conductance. For both 
approaches we have found (rather) closed expressions (valid for a finite channel number) for 
any m-point correlation functions of the transmission eigenvalues, allowing us to calculate the 
first and second moment of the conductance. This complete analytical solution is essentially 
restricted to the unitary case where both Fokker Planck approaches can be solved in terms 
of fermionic one particle Green's functions. At the moment, the extension of this solution 
for arbitrary to the orthogonal and symplectic cases looks difficult, since the Hamilton 
operator obtained after the Sutherland transformation does not reduces to the sum of one 
particle operators. 
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APPENDIX: FOKKER-PLANCK EQUATION FOR THE TRANSMISSION 
EIGENVALUES IN THE 5-BROWNIAN MOTION ENSEMBLE. 

In this appendix, we derive the Fokker-Planck equation ( ^.7[ ) for the case (3 = 2. In order 
to do this, we consider the matrix 



^ V ; ^ V ; 2 yOuiJ \ tsm{2<^) cos(2(^) ^ \^ u\ 



where 



^0 -1 

The matrix Q{t + 6t) then has the form 



(Al) 
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( cos{2{ip + 6ip)) tsm{2{ip + Sip))\ (vl{t + St) \ 

\ism{2{if + 6ip)) cosi2{ip + 6ip)) J [ u\it + 6t) ) ' ^ ' 

On the other hand the evolution equation S(t + 6t) = S{t)e^^-^ imphes 

Q{t + 5t) = S{t)e'^^e-'^^^^^^J:,Sityj:, = S(t)e''^«S,S(t)^S, (A4) 



with 



i6Q = i{6X - S,5XS,) + i \bX, S^(5XS,] + 0[6X^) . (A5) 
The matrix bQ has the form 

where (5m is an arbitrary complex x A^-random matrix with 

{brriij) = , {5rhij5mki) = AD5t 6ikSji , {6mij6mki) = . (A7) 
The matrix S{t) in (S3) contains still the unitary matrices V2 and U2 which can be taken 



into account by the replacement V25mu2 = Sm, This gives for 6m the same statistics than 
those described by Eq. ([A^). 

At this point, we can make an important remark concerning the case where the matrix 
6X obeys a more general distribution than that of Eq. ( |2.6| ). In this case, the property 
( |A7| ) for the matrix 6m is of course no longer valid. But let us assume that the two unitary 
matrices U2 and V2 are independently distributed according to the invariant measure for 
unitary A^ x A^-matrices, i.e. U2 and V2 are described by two independent N x N CUEs. 
Then the average over these two unitary matrices yields just the distribution ( [A7[ ) for the 
transformed matrix 6m = V26mu2- The only information that we need from 6m is therefore 
the invariant quantity 

(tr(W5m)) = ADN^ 6t (A8) 

which relates the time step 6t with the perturbation 6X or 6m respectively. This proves 
that the Fokker-Planck equation derived below is also valid for a more general perturbation 
6X, provided that the two unitary matrices U2 and V2 are for each time t independently 
CUE distributed. The latter assumption is of course of crucial importance. In the case of an 
initial condition of one 2A^ x 2A^ COE or two N x N CUEs for S{0), this assumption is at 
least valid at t = and we assume that it remains correct for arbitrary t > 0. We conclude 
that the Fokker-Planck equation derived for a perturbation 6X obeying ( |2.6| ), still holds if 
one can relate the Fokker-Planck time to the perturbation via ([A^). 

Since we are interested in the statistics of the 6ipj in ( |A3| ) , we have to express the 6ipj in 
terms of 6m. This corresponds to an eigenvalue problem after the transformation 
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Q = Wqu , u 



which imphes 



Q{t) = 


W 


Q{t + 5t) = 


w 






SQ = 


\ 



g2i<^ 

e-2«</' 



\ ,,Q/e^^ 
^ e"*^ J \ e-'^ ^ 

5m) + dm 5m) — 5m \ 
-{5fh) — 5m) —{5m) + 5m) 1 

The matrix W which appears in ( |A10| ) and (|A11| ) is just the unitary matrix 



W = U^ 



vi 

Ml 



U 



(A9) 

(AlO) 
(All) 
(A12) 

(A13) 



where U is defined in ( |A9|) . The matrix W resuhs only in a similarity transformation which 
does not change the eigenvalues. Eq. ( [A^ ) implies that Q(t + 5t) has just the eigenvalues 
^2i{ipj+Sipj) g^j^^ ^-2i(ipj+6ipj) ^ which can be obtained by second order perturbation theory 



VjkVkj 



E 



yj,k+Nyk+N,j 



(A14) 



The perturbing matrix V in Eq. ( A14 ) is just 



V 







\ 

Eq. ( |A14|) leads after some algebra to 

1 
2 



e'^ 
e-'^ 



(A15) 



5ipj 



(A16) 



We insert Eq. (|A12D in Eq. (|A16|) and apply the average ([A7D for the matrix 5m that results 
in 



with 



fj{ip) = 2 J2 {cot{ipj - ipk) + cot(v?j + (fik)) +2cot{2ipj) 

The corresponding Fokker-Planck equation has the form 

( d 



dp{(p, t) _ D ^ d 



dt 



4 ^ difj \\dipj 



(A17) 
(A18) 

(A19) 



which is just Eq. ( p.7| ) for /? = 2 because fj{^) = In -^2(^9) with -^2(^9) given by Eq. 
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FIGURE CAPTIONS 



Fig. 1 Two chaotic ballistic cavities coupled with a n channel contact and each of them 
connected to an electron reservoir through a channel contact. 

Fig. 2 Density of the variable u = (arsinh(A))^/^ for the crossover 2xCUE CUE with 
= 1 at different values of t/tc = 10.0, 1.0, 0.1, 0.01. The curves are calculated 
from Eqs. {^J§ and (lOTl) . 

Fig. 3 Density of the variable u = (arsinh(A))^/^ for the crossover 2xCUE CUE with 
= 5 at different values of t/tc = 10.0, 1.0, 0.1, 0.01. The curves are calculated from 
Eqs. (|3.19|) and (|3.21| ). In addition, for t/tc = 0.1, 0.01 the Laguerre approximation 
is included. 

Fig. 4 Density of the variable u = (arsinh(A))^/^ for the disordered line at three different 
length scales, L/^ = 0.2, 1.0, 10.0. The curves are obtained from a numerical evalua- 
tion of Eq. (PI). 

Fig. 5a Same as Fig. 4 for L/C, = 0.2 but with a modified scale. 

Fig. 5b Same as Fig. 4 for L/C, = 1.0 but with a modified scale. 

Fig. 5c Same as Fig. 4 for L/C, = 10.0 but with a modified scale. The Gaussian approximation 
(Eq. ( [4.38| ), dotted line) is included for comparison. 
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